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1. Introduction 



In a quite recent work Aspelmeier and Moore [Q have considered the sample-to-sample 
free energy fluctuations in finite dimensional spin glasses via the replica method. To that 
effect they reconsider higher order terms in the replica number n and they conclude that 
the Parisi symmetry breaking scheme || does not give the correct answer for these higher 
order terms. Finally they propose a modified symmetry breaking scheme that resolves the 
problem. 

What we set out to do here is the following. Starting from the same truncated Hamil- 
tonian (AM. 3) we solve exactly the discrete stationarity equations for R steps of replica 
symmetry breaking, namely we obtain the R + 1 values of q^p indexed by their overlap 
values qo,q±, •••,<?/? (together with q aa = qR+i = 0) and the R values of Parisi box sizes 
Pl,p%, ■■■iPr together with the two fixed boundary values po = n and -Pr+i = 1. As a 
result, we find two families (a), (b) of solutions associated with two possible values of qo, 
namely, letting g = w/(2y) 

(a) q = ^g 

In this case the corresponding free energy is identical to the Kondor || result 

Qr?f> 

n/<°>(n)=n/-— wg 3 (1.1) 
o40 

(b) q = 

The free energy is now larger 

nf( b \n) = nf (1.2) 

The (b) solution is therefore the appropriate one to choose, both solutions having a 
non-negative Hessian spectrum when R — > oo. Among the family of solutions (b) with 
qo = and free energy f^ b \ we will pick a reference solution with a set of values qt,Pt, 
t = 1,2, ...,R. All the other solutions will be shown elsewhere Q to correspond to a 
(discrete) reparametrization for large R. With that set of values, we proceed and compute 
the contribution to fluctuations, with a result that matches for R —>■ oo the Aspelmeier 
and Moore ones Q. We thereby establish that there is indeed no need for modifying the 
Parisi replica symmetry breaking scheme. 
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2. Solution of the stationarity equations 



The stationarity equations are derived from the free energy functional 

t=Q 

where we have used the replica Fourier transform q of q 

R+l R 

Qk = Y^ Pt ( qt ~ qt ~^ = ^2p*(lt - Qt-i) - QR (2-2) 

t=k t=k 

Combining the stationarity equations, we obtain in the end 

9Pt = \{qt + qt-i) t=l,2,...,R 

(q t - q t -i) 2 = (qt-i ~ qt-2) 2 = ■■■ = (<?i ~ qo) 2 t = l,2,...,R 
Here we concentrate on the particular reference solution such that 

q t - q t -i = q t -i - q t -2 = •■• = 91 - <7o = — p — (2.4) 

which leads to 

t 

qt = qo + {qR-qo)^ t = 0,l,...,R 



, 2t - 1 

9Pt = qo + {qR - Qo) R t=l,2,...,R 



(2.5) 



together with qR+% = 0, Pr+i = 1. Besides one has two more equations that determine qo 
and q R 

E(, R )- i(^) 2 = (2-6) 
where E(qn) = r — wqR + uq R , and which is valid for R > 0, and 

go (^(g«) + f (3<7Po - 2<fo)) = (2.7) 



valid for all i?. Note that if R = 0, g# = qo, then ( |2.5|) is a tautology and only ( |2T7| ) 
survives, leading to the standard result wq = 2r/ (2 — p ) +0(r 2 ). In fact, one is interested 
in the limit of large R, whereby (|2.6|) yields the relationship 

E{q R ) = (2.8) 

and from ( |2.7[ ) either q = or q = 3gpo/2 as respectively in the cases (b) and (a). Note 
that q t is monotonous except for its last step (qR+i = 0), and p t is monotonous except for 
its first step (when p is kept fixed at a value n^0). 

In the continuum limit, where t/R — > x and — > </(x), — >■ p(x), we get for x in the 
open interval (0, 1) 

q(x) = gp(x) = gqRX < x < 1 (2.9) 
We now proceed to get the fluctuation contribution as in (AM. 5). 
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3. Fluctuations: the Replicon sector 

We have as in (AM. 9) 

nSf Rep = ^ J ^-^I Rep (p) (3.1) 

where 

R R+l 

I Rep (p)=nJ2 Kr;kJ)Log(p 2 + \(r;kJ)) (3.2) 

r=0 k,l=r+l 

Here the Replicon eigenvalue A is 

X(r;k,l) = -2r - wq k - wqi -2yq 2 (3.3) 
The multiplicity fi(r; k, I) is given by 

M r ; fc > = \{Pr- p r +l)fJ.(k)(i(l) (3.4) 

Wh6re f^-^ k>r + l 

m = \ Pk ^~ 1 k = r + 1 (3-5) 

We note that po is absent from q, since even if the index k, I were allowed to take the value 
0, it would appear in the vanishing combination poqo- The po dependence can therefore 
only arise from the multiplicity. Collecting the terms in po we get 

R+l R R 

T U 

J-Rep — 



^M^(OLog(p 2 + A(0;fc,/)) + ^p r £ ^( ^tcr-l^M) 



(3.6) 



(3.7) 



With the use of (|2~5|)-(p^D and we get 

A(r; fc, /) = -2E(g B ) + 9^)\\{k 2 + I 2 ) - r 2 - (k + I) + 1) 

We note that the lowest (Replicon) eigenvalue is given by 

A ( r; r + l,r + l) = -i(^) 2 r = 0, 1, 2, R (3.8) 

hence we find an instability, except in the limit R — > oo where this is suppressed. All other 
eigenvalues are positive. We thus have in the Parisi limit R+l zero modes arising from 



the negative eigenvalues (3^8). In that limit, one has 



I 



Rep 



\ j 1 xdxd x { jf * fd k £ fd^p 2 + gf(k 2 + l 2 - 2x 2 )) } 



which coincides with (AM. 12). 



(3.9) 
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4. Fluctuations: the longitudinal-anomalous (LA) sector 

We now have 



V f d D p 

R+l 



nSf LA = — I tttZ\d SI la(p) 



SIla(p) = n^^(/c)LogdetA fc (r,s) (4.1) 



fc=0 



A fc (r, g )=e- W ^ (r " ) ^" 1) 
where 5 Kr denotes the Kronecker delta, while we have 



„(*-!) 

Irs 



Ps ' -Ps+1 

Ps S > k — 1 

2p s s < jfc - 1 (4.2) 



. ( r ) _ f p 2 + A(r; r + 1, k) k > r + 1 
fc 1 p 2 + A(r; r + 1, r + 1) fc < r + 1 



Expanding the determinant A& (r, s) yields 



detA*(r, S )=i+x;(-^) n e n^-^-jw- ( 43 ) 

m=l 0<si<s 2 <...<s m i=l 

where we have set so = 0. In order to have po occurring in the determinant, i.e. in one 
of the 5 S ^ 1 , we need Sj = hence the only possible term is si = 0, but the prefactor 
q Sl — Qs = Q Sl = Qo vanishes and there is no po contribution form the Log again. The 
only po contribution comes from the multiplicity np(k) which now cannot sustain an n 2 
contribution!!. 



5. Conclusion 

With no contribution to fluctuations from the LA sector, we conclude that the full 
answer is given by the contribution from the Replicon sector as of ( |Q| ) and (|3.9|) , thus 
corroborating the result of M. 



3 One may then ask what becomes of the term k = which has a factor n/p<y. It is actually 
given by SI la — L°g det(Ao(r, s)/Ai(r, s)), but with 5^ = Si ^ —poS^o, whereas Ai (s) = Ao(s). 
Again, SI^a reduces to a contribution ~ poS^o which vanishes with go- 
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6. Extension to Aspelmeier, Moore, Young calculation 

There the authors give an analytic answer to a long standing problem: computing 
the interface free energy of the Ising spin glass. They show that (part of) the answer is 
obtained by computing Z n Z m and, in the associated free energy, leaving aside the terms 
in (n + m), (n + m) 2 and keeping the terms in nm. From ( |3.6j ) we clearly see that, leaving 
aside a (n + m) 2 contribution, one is left with a term 

Ipip) = IRepip) = -nm M*OMOLog(p 2 + A(0; k, /)) (6.1) 

k,l>l 

replacing the quadratic term of ( |3.6| ) and associated with periodic boundary conditions. 
The fluctuation contributions from the off-diagonal blocks of the Hessian (mixed sector, 
associated with antiperiodic boundary conditions) is much easier to deal with since it is 
exactly given by 

R+l 

Iap(p) = +nm ^ u(k) fi(l)Log(p 2 + A(0; k, I)) (6.2) 

kd=0 

eigenvalues and multiplicities matching the R = calculation presented earlier by the same 
authors 0. The LA sector is represented here respectively by k = I = and k = 0, I > 1 
or k > 1, I = 0, the Replicon sector by /c, I > 1 as in (|BTT|) . But with qi = qo we now see 
that the two contributions periodic ( |6.1|) and antiperiodic ( |6.2|) are formally identical but 
for their sign. This result again matches the one proposed in ||. 
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